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Abstract The inverse solution procedure that enables the 
identification of the defect position in a beam from the reso¬ 
nance frequency was exploited. Resonance frequency shifts 
of a power spectrum due to defects in a longitudinally 
vibrating beam when both ends are free were investigated 
by both numerical and experimental analysis. Calculation 
by a transfer matrix method showed that the frequency shift 
was large when the defect position coincides with a node of 
vibration and that no shift occurs when it coincides with a 
loop of vibration. The frequency shift could be approxi¬ 
mated by a sinusoidal curve. Calculation results agreed well 
with those of the experiment in which artificial round holes 
were drilled as the defect model. Experimental equations 
predicting the amount of the frequency shift in function of 
the defect position were obtained. In the inverse procedure, 
the defect position was determined by comparing the reso¬ 
nance frequencies between the experimental and estimated 
power spectra so that the coincidence factor S(x) became a 
minimum. The results showed the validity of the proposed 
method to identify the defect positions of fewer than two 
predominant defects. 

Key words Longitudinal vibration • Power spectrum ■ Iden¬ 
tification of defect • Transfer matrix method • Inverse 
problem 


Introduction 

Existence of a defect that affects the strength of a structural 
timber such as a knot induces a resonance frequency shift 
of vibration. 1-3 Numerical calculation using a finite-element 
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method has been applied to the problem. 3 Several ideas that 
aim for the identification of a local defect have been pro¬ 
posed, for example, the flexural curve of a beam in static 
bending test 4 and the flexural curve associated with a trans¬ 
fer function in flexural vibration. 5-7 In the Japanese wood 
industry, the major type of stress grading machine for sawn 
timbers for Japanese wooden construction has recently 
shifted into the longitudinal vibration type machines using 
a tapping excitation, because the dimension of the struc¬ 
tural timber in traditional Japanese construction is larger 
than that of the 2 by 4 foot timbers used in wood construc¬ 
tion in North America. 

The merit using the power spectrum of the vibration is 
that the spectrum can be obtained from a transient sound 
caused by one hit of a hammer. This method is easily intro¬ 
duced into a current stress grading machine of the longitu¬ 
dinal vibration type because the frequency data have been 
supplied from a fast Fourier transform (FFT) analyzer 
included in the dynamic stress grading machine. In this 
article, the effects of defects on the frequency shift of the 
power spectrum were estimated by means of a transfer 
matrix method. Moreover, experimental work was con¬ 
ducted by introducing artificial round holes as structural 
defects. Then, the prediction of defect position was exam¬ 
ined. The defect position was inversely identified from the 
coincidence test of the power spectrum between the experi¬ 
mental data and calculation by the transfer matrix 
method. 


Theory 

Transfer matrix method 

The transfer matrix method is a practical tool for solving 
vibration problems of nonuniform cross-sectional bodies. 8 
This tool has been used in the vibration analyses of wood 
members such as slender logs, 9 plywood-webbed box 
beams, 10 the sound boards of a marimba, 11 and the flattening 
of a resonance frequency shift by defect. 9 In the transfer 
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matrix method, 8 a body is considered to be a linear combi¬ 
nation of spring elements and point mass elements. Under 
the boundary condition of both ends free, the solution with 
respect to frequency is the resonance frequency to be 
obtained. In the calculation, the effect of the defect is evalu¬ 
ated as a decrease of a Young’s modulus of the correspond¬ 
ing local field matrix. Here, Hatter’s FORTRAN program 
of the torsional vibration problem 8 was rewritten into a 
BASIC program and modified to adopt it for longitudinal 
vibration of a beam subjected to the both ends free condi¬ 
tion. A similar numerical calculation using a finite-element 
method has been applied to the same problem of frequency 
shift by defects. 3 


Identification of defect position 
General strategy 


When the defect position is known, the resonance 
frequencies are given from the calculation results, as shown 
later. However, identification of the defect position 
from the power spectrum is a complicated problem. Our 
strategy is based on the following method of least squares. 
The defect position x is determined so that the differences 
of the resonance frequencies between the experimental 
and estimated spectra become minimum. The following 
coincidence factor, S(x), is defined as a parameter of 
identification. 

A relative frequency shift of an arbitrary vibration mode, 
rsf n , due to a defect is given by 


rsf n = - lr " x 100 = fr " nlr ' a x 100 (%) 

frn, 0 n frl,Q 


( 1 ) 


where fr n is an experimental resonance frequency of an nth 
mode of a specimen with a defect, and then fr n 0 and fr l0 are 
resonance frequencies of an nth and the fundamental modes 
in a defect-free condition, respectively. In our calculation 
model, the following relative frequency shift k n (x) is given 
as a function of the defect position x(0 < x < 1): 


k n (x) = 


K{x)-nfr[ Q 
nfrlo 


x 100 


( 2 ) 


where fr'Jx) and fr j 0 are the resonance frequencies of the 
nth mode given from the result of the transfer matrix 
method later. The value x that gives the minimum value of 
5(x) in Eq. 3, the coincidence factor, is the optimal solution 
of the defect position. 

S (x) = {rsf\ - /q(x)} 2 + {rsf 2 - k 2 (x)} 2 + • ■ • + 

{rsf N - Mx)} 2 = - k n(x)} 2 (3) 

n =1 


where N is the highest vibration mode to be considered. In 
this article, N= 5. A piece of information given by one reso¬ 
nance frequency may be too poor to identify the defect 
position, but integration of information of several vibration 
modes would increase the degree of confidence on the iden¬ 
tification of the defect positions. 


Coincidence factor in a two-holes case 

Assuming an additional low in the frequency shift, the 
whole frequency shift of the nth vibration mode Rsf n when 
the number of holes is two is given by the summation of the 
effect of each hole. 

Rsf n =rsf^ + rsfl 2) (4) 

where the superscripts (1) and (2) are the number of the 
two holes. 

The verification of this assumption is proved in this 
experiment later. 

The coincidence factor in the two-holes system is 
extended from Eq. 3 as follows: 

5(Xi, X 2 ) = X [ Rs fn - {K (*1 ) + kn (x 2 )}] 2 (5) 

n =1 

where k n (xf) and k n (x 2 ) are the calculation values of relative 
frequency shifts in the nth vibration mode to be identified 
and x, and x 2 are the positions of the two holes. In a speci¬ 
men having two holes, the following identification in two 
steps was applied. At the first step, the first defect position 
x' is determined so that the coincidence factor becomes 
minimum, supposing a single-hole system, and then k n (x\) 
is given by Eq. 2. Then, the second hole is searched by cal¬ 
culating the following coincidence factor S(x 2 , x 2 ) given by 

S (x, ,x 2 ) = X [{Rsfn - k n (x?)}- k n (x 2 )f (6) 

«=i 

For more than three holes, Eqs. 4 and 5 can be extended to 
a general case by the same procedure. 


Materials and methods 

The effect of defects on the resonance frequency shift in the 
transfer matrix method was calculated for the condition that 
the degree of reduction of Young’s modulus at the field 
matrix of the defect positions was 50% from the defect-free 
Young’s modulus. Here, the reduction is equivalent to a 
cross-sectional loss. The number of the field matrices of the 
beam was chosen to be 67 to adopt the defect size of the 
following experiment. Then, the field matrix was replaced 
by the modified element successively from the end element. 
The following material constants were used: Young’s 
modulus, 10 GPa; density, 0.5 g/cm 3 ; length, 1 m; and cross 
section, 4.5 cm 2 . 

The number of elements greatly affects the results. To 
check the accuracy of calculation, the resonance frequency 
of a homogeneous beam calculated by the transfer matrix 
method was compared with its analytical solution of vibra¬ 
tion theory for cases of different numbers of field matrices 
from 3 to 150. The calculation error of the resonance fre¬ 
quency was less than 0.001% for 67 field matrices. 

Eight spruce beams (Picea sp.) with dimensions of 
1000 x 30 x 15 mm (longitudinal, tangential, and radial 
directions) were used for experimental verification. The 
specimen density and the oven-dry moisture content were 
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0.44 g/cm 3 and 8.8% on average, respectively. Round holes 
were drilled through the wide faces of each specimen. The 
diameter of the hole was 15 mm, that is, half the width of 
the specimen. Four holes were drilled finally at 1/8, 1/4,1/3, 
and 1/2 positions of the specimen length, but the order of 
drilling was different for the four pairs of specimens. The 
data of each pair were averaged in the following analysis. 
Additional verifications were conducted for a piece of com¬ 
mercial structural 2 by 4 foot dimension lumber with knots. 
This lumber included a large knot whose knot ratio was 
around 0.3 at the relative distance 0.45 from an end of the 
lumber. 

A FFT spectrum analyzer was used to obtain power 
spectra. Specimens were supported on pieces of foamed 
polyurethane. FFT utilized the standard function of the 
Cooley-Tukey type algorithm of the FFT spectrum ana¬ 
lyzer. The frequency resolution increases with increase of 
the number of vibration modes. The maximum error in the 
fundamental vibration is estimated to be around 1% and 
those of the overtones are smaller. To diminish the restric¬ 
tion of free vibration by the support system, the odd-num¬ 
bered vibration was measured by supporting the specimen 
at the center and the other vibration mode was measured 
by supporting it at the nodal positions of the second-num¬ 
bered vibration mode. 


Results and discussion 

Simulation by the transfer matrix method 

Figure 1 shows the resonance frequency shifts from the 
fundamental to the fifth modes. The calculated values were 
symmetrical with respect to the center position. Then, only 
the values over the half-span of the beam between 0 < x < 
0.5 were shown. In each vibration mode, the frequency shift 
is large when the defect position coincides with the node of 
vibration, and no shift occurs when it coincides with the 
loop of vibration. The magnitude of the maximum fre¬ 
quency shift was 1.47%. 


0 0.1 0.2 0.3 0.4 0.5 



Fig. 1. Relative shift of resonance frequency by a hole calculated from 
the transfer matrix method. Solid lines show approximation by sinusoi¬ 
dal curves in Eq. 8. Filled circle. First mode; open circle, second mode; 
filled triangle, third mode; open triangle, fourth mode; cross, fifth 
mode 


The frequency shift at an arbitral defect position k n (x) 
was approximated by a sinusoidal curve. 

kJx)= ( ^cos{2nnx)-^ (7) 

where a is a maximum amplitude of the frequency shift and 
n is the number of vibration mode. 

The approximated curves by cosine functions are also 
shown by solid lines in Fig. 1. The curves give a good fit to 
the calculation results by the transfer matrix method. The 
pooled data gave the following equation: 

k n (x) = ^^-{cos(2nnx)-l}. (8) 

Effect of the positions and number of holes on the 
frequency shift 

Figure 2 shows the effect of hole position on the frequency 
shift of the spectrum peak when four holes are drilled suc¬ 
cessively in the order of the 1/2, 1/3, 1/8 and 1/4 positions 
of the beam length. In the first drilling at the 1/2 position, 
peak shifts of the odd-numbered vibration modes occurred. 
However, those of the even-numbered modes did not occur. 
In the second drilling at the 1/3 position, additional fre¬ 
quency shifts occurred at the first, second, fourth, and fifth 
numbered modes. However, the shift of the third numbered 
vibration mode did not occur. In the third drilling at the 1/8 
position, the peak shift of the fourth numbered vibration 
mode occurred greatly, and large shifts occurred at the 
second and the third numbered vibration modes; that of the 
fifth numbered vibration mode was slight. In the last drilling 
at the 1/4 position, a large shift occurred at the second 
numbered vibration mode, and slight shifts occurred at the 
odd-numbered vibration modes. However, no shift occurred 
at the fourth numbered vibration mode. These results show 
that a large frequency shift occurs when the hole position 



Fig. 2. Effect of hole position on relative frequency shift. Four holes 
are drilled successively in the order of 1/2,1/3,1/8, and 1/4 positions of 
the beam length from an end 
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Fig. 3. Relationship between relative frequency shift by the experi¬ 
mental rs/ stogle and that by the calculation k n ( x) 


coincided with the node of each vibration mode and that 
the shift did not occur when it coincided with the loop of 
each vibration mode. In the other cases, frequency shifts 
occur more or less according to the position of defects. A 
slight increase of the resonance frequency is shown at the 
1/2 position of the even vibration modes. A loss of mass at 
a hole position can induce an increase of the resonance 
frequency for a combination of the hole position and the 
vibration mode. The same phenomenon has been reported 
in flexural vibration by Nakayama. 1 

Prediction of relative frequency shifts 
Single-hole case 

Figure 3 shows the relationship between the relative fre¬ 
quency shift by the experimental rsf siTtglt and the calculated 
k n (x). The following regression line was obtained for pooled 
data from the fundamental to fifth ordered vibration modes 
of single-hole specimens: 

rs/; ingle = 3.173k„(x) (R 2 = 0.89) (9) 

The correlation coefficient of determination, R 2 , was 
very high, but the magnitude of the calculation value was 
about 30% of the experimental value rs/ sillgle ; this may be the 
result of the effect of stress concentration around a hole, 
which was ignored in the transfer matrix calculation. The 
same effect has been reported in flexural vibration by 
Nakayama. 1 

Multiple-hole case 

The following test was done to verify the additional low on 
the relative frequency shift for multiple holes. The estima¬ 
tion of the relative frequency shift of the single hole rs/ single 
is given by substituting Eq. 8 in Eq. 9: 

^/single = 2.332(cos2«7tx -1) (10) 



Estimated relative frequency shift (%) 

Fig. 4. Relationships between relative frequency shift by experimental 
Rsf mu in and by the additional low based on the data of a single hole 

,cal 



Relative distance from an end 

Fig. 5. Identification of a single hole drilled at 1/2,1/4, and 1/8 positions 
of the beams 


The relationship between the relative frequency shift by the 
experimental Rsf multi and by the additional low based on the 
data of the single-hole Rsf mmjai i is shown in Fig. 4. The data 
were pooled from those of the first to fourth holes and the 
fundamental to fifth vibration modes. 

Rsf mm = 0.8504 Rsf multiM , R 2 = 0.767 (11) 

The experimental values are 85% of the calculation values, 
applying the additional low, and the coefficient of determi¬ 
nation is sufficiently large. The linearity in Eq. 11 proves 
the validity of the additional low in the frequency shift. 

Identification of the defect position 
Application to the single-hole case 

Figure 5 shows the results when the single hole was drilled 
at 1/2, 1/4, and 1/8 positions of the beams, respectively. 
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Relative distance from an end 


Fig. 6. Identification of two holes drilled at 1/2 and 1/3 positions of the 
beam 



0 0.1 0.2 0.3 0.4 0.5 

Relative distance from an end 

Fig. 7. Identification of two holes drilled at 1/2 and 1/8 positions of the 
beam 

Minimum values of the coincidence factor S(x) were located 
at each hole position. The coincidence factor was suffi¬ 
ciently small, and clear identification of the hole position 
was shown. 

Application to the two-holes case 

Figures 6 and 7 show the results when two holes were drilled 
at 1/2 and 1/3 positions and 1/2 and 1/8 positions of the 
beams, respectively. These holes were identified by two 
steps. In Fig. 6, the first hole was detected at the 1/2 position 
because the confidence factor was minimum at the center 
of the beam. Applying Eq. 6, the second hole was detected 
near the 1/3 position of the beam. In Fig. 7, the first hole is 
detected at the center and the second one is searched at the 
0.17 position of the beam. The second hole was searched as 
slightly inside the 1/8 position. The multiple-step identifica¬ 
tion worked efficiently. However, the resolutions of the 
third holes were not sufficient, because the coincidence 
factor became rather larger than those of the second 
searches. Thus, the number of holes to be identified was 
limited to two. 



0.00 0.10 0.20 0.30 0.40 0.50 

Relative distance from ends 

Fig. 8. Location of knots from ends in the 2 x 4 in. lumber through the 
length 

Application to structural lumber having knots 

Commercial lumber may include defects such as knots. The 
resonance frequencies of all vibration modes can be shifted 
by existence of knots. Therefore, a difficult problem in this 
case is “What is the fundamental resonance frequency in 
the defect-free condition of each piece of lumber?” Here, 
the following maximum frequency in the normalized reso¬ 
nance frequency frjn among the five vibration modes was 
chosen as a standard frequency, fr*. 

/r*=max(/r„/«; n = 1,2,3,..., 5) (14) 

rsf: = frJn f ~ fr * x 100 (15) 

This procedure was applied to a piece of 2 x 4 in. struc¬ 
tural lumber that has a large knot near the center, whose 
knot ratio was 0.31 at the relative distance 0.45 from an end 
of the lumber. Figure 8 shows the location of knots through 
the length in the same piece of lumber. The relative fre¬ 
quency shifts of the standard frequency at the first to fifth 
vibration modes were -3.147%, 0%, -3.263%, -1.224%, 
and -2.098%, respectively. The frequency shifts of the first 
and third modes were large, that of the fifth mode was 
moderately large, and those of the even modes were small. 
This result means that a predominant knot will locate near 
the center of the lumber. In Fig. 9, the knot was searched 
at the relative distance 0.45. In the second search, the coin¬ 
cidence factor became rather larger than in the first search, 
which means that the second defect is neither predominant 
nor negligible compared with the first defect. 


Conclusions 

The inverse solution procedure that enables the identifica¬ 
tion of the defect position from the power spectrum of the 
longitudinal vibration was exploited. It is composed of the 
estimation of the frequency shift, which was calculated by 
the transfer matrix method, and the coincidence test, which 
searched the optimal solution between the experiment and 
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Relative distance from an end 

Fig. 9. Identification of knot location from an end in the 2 x 4 in. 
structural lumber 


the calculation. The defect position was well identified, but 
the procedure was limited to the identification of two domi¬ 
nant defects. 
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